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e + e~ — > T7 + 7r~e + e~ : a potential background 
for cr(e + e~ — > ir + 7r~) measurement 
via radiative return method * 
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A Monte Carlo generator (EKHARA) has been constructed to simulate 
the reaction e + e~ —> n + Tr~e + e~ based on initial and final state emission of 
a e + e~ pair from e + e~ — > 7r + 7r _ production diagram. A detailed study of 
the process, as a potential background for er(e + e~ — » -k + ii~) measurement 
via radiative return method, is presented for <£>- and B- factory energies. 

PACS numbers: 13.40.Ks,13.66.Bc 

1. Introduction 

The radiative return method |.P (look j2J for a short introduction) is 
a powerful tool in the measurement of a(e + e~ — > hadrons), crucial for 
predictions of the hadronic contributions to a^, the anomalous magnetic 
moment of the muon, and to the running of the electromagnetic coupling 
from its value at low energy up to Mz (for recent reviews look [31 HJ [3] ) . 
Due to a complicated experimental setup, the use of Monte Carlo (MC) 
event generators ^ H3 |S1 , which includes various radiative corrections 
jlUl lllj is indispensable. Some more extensive analysis of that subject can 
be found also in |12l 1131 I14j . The most important hadronic mode, i.e. 
7r + 7r~, is currently measured by KLOE ^3^3 El dU> and BaBar ^5] by 
means of radiative return method. This measurement can suffer from a 
background from the process e + e~ — > 7r + 7r _ e + e _ , as suggested in J2Uj, for 
at least two reasons: 1. At present KLOE measures only pions (+ missing 
momenta) in the final state and for that particular measurement there is 
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no difference between photon(s) and pair production. 2. The e + e~ pair 
can escape detection, being lost outside a detector, e.g. in the beam pipe, 
or having energy below a detection threshold. Again a Monte Carlo study 
is unavoidable, if one likes to know the actual value of the pair production 
contribution in a given experimental setup, as the analytical, completely 
inclusive, calculations might lead to misleading results. 



2. Monte Carlo program EKHARA and its tests 




Fig. 1. Diagrams contributing to the process e + {p\)e~(p2) — ► 

7r + (7ri)7r~(7T2)e + (qi)e~(<72): initial state electron pair emission (a), final 
state electron pair emission (b), pion pair emission from t-channel Bhabha process 
(c) and 7*7* pion pair production (d). 

In Fig^ different types of diagrams contributing to process e + e~ — > 
7r + 7r~e + e~ are shown schematically. In the present version of the Monte 
Carlo program we include only two gauge invariant sets of diagrams from 
Fig^i and^>. The former represents initial state radiation (ISR), and the 
latter final state radiation (FSR), of an e + e~ pair from e + e~ — > ir + ir~ 
production diagram. We use scalar QED (sQED) to model FSR e + e~ pair 
emission and p dominance model for 7*(p*)7T7r coupling (see for details). 
The diagrams from Fig^, representing pion pair emission from t-channel 
Bhabha process, together with s-channel Bhabha pion pair emission (not 
shown in FigQ), will be included in the new version S2li of the presented 
generator, completing the discussion of this paper. The contribution from 
7*7* pion pair production process (Fig^i) is negligible for DA<1>NE energy 
|17| . and as its interference with other diagrams does not contribute to the 
cross section integrated over charge symmetric cuts, these contributions are 
not relevant, at least for ^-factories. 

For parametrisation of the phase space we use the following variables: 
Q 2 = (tti + 7T2) 2 -invariant mass of -K + ir~ system, k\ = {q\ + g2) 2_ i nv a r iant 
mass of e + e~ system, polar and azimuthal angles of Q momentum, defined 
in the centre-of-mass (CM) frame of initial e + e~ pair, polar and azimuthal 
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angles of tti momentum, denned in Q- rest frame and polar and azimuthal 
angles of q\ momentum, defined in k\- rest frame. All four vectors are 
boosted into the initial e + e~ CM frame after being generated and all nec- 
essary cuts can be applied at this stage of the generation. Multi-channel 
variance reduction method was used to improve efficiency of the generator 
and all details will be given in a separate publication |21j . 

We have performed a number of 'standard' tests to ensure that the 
written FORTRAN code is correct. Gauge invariance of the sum of the 
amplitudes was checked analytically separately for set of diagrams from 
Fig^i and^>. We use helicity amplitudes in EKHARA to calculate square 
of the matrix element describing the e + e~ — > 7r + 7r - e + e~ process, but as a 
cross check, we have used also the standard trace technique to calculate the 
square of the matrix element, summed over polarisations of initial and final 
leptons. Both results were compared numerically scanning the physical 
phase space, and the biggest relative difference between the two results 
found was at the level of 1CT 9 . It was necessary to use quadrupole precision 
of the real numbers for the trace technique result, as one can observe severe 
cancellations between various terms. The phase space volume calculated by 
Monte Carlo program was cross checked with the Gauss integration and the 
relative difference at the level of 1CP 5 was well within the errors of the MC 
result, which were of the same order. 
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Fig. 2. EKHARA results compared with analytical results of [221 ( see text for 
details). The errors come from MC integration. 

Inclusive analytical formulae from [22] provide additional, nontrivial 
tests of the implementation of the contributions from Fig^i. Formula (23) 
from [22j provides Q 2 differential cross section (other variables are inte- 
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grated over the whole allowed range) valid for large Q 2 . In FigEKj we 
compare the values of the integrals, over 10 equally spaced intervals of Q 2 , 
obtained by means of MC program and one-dimensional Gauss integration 
of the above mentioned analytical formula. The Gauss routine, which we 
use, guarantees precision of 12 significant digits. One observes a relatively 
good agreement for values of Q 2 ~ s and worse for Q 2 nearby tt + it~ produc- 
tion threshold, as expected from the applicability of the analytical formula. 
EKHARA results agree much better (see FigEJo) with known analytically 
doubly differential cross section d( ^ dk z |H3 HU , integrated over the whole 

allowed range of k 2 and 10 equally spaced intervals of Q 2 . The exact ana- 
lytical result was integrated numerically, using recursively one-dimensional 
8-point Gauss procedure and dividing the region of integration into pieces 
small enough to guarantee the overall accuracy of 10 significant digits. From 
FigEJa it is clear that a technical precision of EKHARA of the order of 0.1% 
was achieved. 



3. Monte Carlo Data Analysis 





a) 


75 = 1 
no eut 


.02 Ce\ 


I 4 




s 










* 
























* 












+ 

+ 












* 
+ 










* 


* 

* 








* 
* 
* 










f 







0.2 0.4 0.6 O.f 

Q 2 (GeV 2 ) 




0.4 0.6 0.8 

Cf (GeV 2 ) 



Fig. 3. Comparisons of ISR and FSR contributions to e + e — > ir + ir e + e cross 
section at DA$NE energy. 

For ISR of e + e~ pair ( Fig^i), a factorisation similar to photon emission 
holds [22] and adding ISR pair production to ISR photon production results 
just in a change of the radiator function, thus radiative return method 
still can be used 0. On the other hand, FSR of e + e~ pair ( FigOa) is 
model dependent, the same way as it is the emission of a real photon, and 
the question of its relative, to ISR, contribution to the cross section is as 
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important as for the photon emission. One can observe, that e + e~ pair 
emission resembles a lot photon emission, with big contributions of FSR 
for inclusive configurations (Fig|3K) of a <£-factory, which can be easily 
reduced, by suitable cuts, to a negligible level (FigEJa). Moreover, the cuts 
which reduce photon FSR reduces as well the e + e~ pair FSR. In addition, 
analogously to photon FSR, e + e~ pair FSR is completely negligible at B- 
factories. 
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Fig. 4. Ratio of differential cross sections of the process e + e — » 7r + 7r e + e (pair) 
and e + e~ — » tt + tt~ + photon(s) (ph). 

The most relevant information, how big is the contribution of the process 
e + e~ — > 7r + 7r~e + e~ in comparison to the main process used in the radiative 
return method, mainly e + e~ — > 7r + 7r~ +photon(s), is presented in Fig0]for 
DA<l?NE energy, both without any cuts (Fig^t), and with cuts resembling 
KLOE event selection QU Q2] (Figgja). The results of e+e" -> vr+vr" + 
photon(s) cross section were obtained using PHOKHARA 3.0 MC generator 
[Hj and in the following, whenever we refer to e + e~ — > 7r + 7r~ + photon(s) 
cross section we mean cross section obtained using PHOKHARA 3.0. As 
one can see from Fig^J the contribution of the e + e~ pair production is 
below 0.5%, independently on the cuts. It is Q 2 dependant, being big 
for low Q 2 values. Even if it is small, this 0.5% contribution can become 
important, when aiming at the precision below, or of the order of 1%, for 
the e + e~ — > tt + tt~ cross section measurement. At S-factories, the relative 
contribution of the pair production might be as big as 0.9% (FigJS^) and it 
is again Q 2 and cut dependent (FigEh)- 

As stated already, ISR of electron pairs can be treated in a similar 
way as ISR of photons resulting in the change of radiator function in the 
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Fig. 5. Ratio of differential cross sections of the process e + e — > tt + 7t e + e (pair) 
and e + e~ — » 7r + 7r~ + photon(s) (ph). 



radiative return method. However, one can alternatively try to treat it as a 
background to the process with photon(s) emission. In this case, the most 
natural way of reducing that background is to veto the electron (positron) in 
the final state. In Fig|f)]we show an example of such a procedure performed 
for ^-factory energy. We assume here that an electron or positron can be 
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Fig. 6. Non-reducible pair production background at DA$NE energy : (a) no 
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< 15° or 9 w + +7r - > 165° and M tr cut. For (a) and 
.+ < 20° or 9 e + > 160°) and (0 e - < 20° or 



50° < 6> ± < 130°, 6 



7T++7T 



(b) c2 denotes additional cuts: 

B e - > 160°). 



ENowakUstronOS'new printed on February 1, 2008 



7 



seen, and the event rejected, if its angle with respect to the beam axis is 
bigger then 20°. FigUJi shows that up to 50% of the events pass the rejection 
procedure, when no other cuts are applied. However, in the case of KLOE 
event selection, which requires that the e + e~ pair is emitted along the beam 
axis, one rejects only a small fraction of these events (FigEb). 
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Fig. 7. Non-reducible pair production background at _B-factory energy : (a) no 
cuts in (^p)tot; (b) for both (^?) C 2 and {-^)tot cuts on pion angles are imposed: 
50° < 0^ < 130°. For (a) and (b) c2 denotes additional cuts: (0 e + < 20° or 
9 e + > 160°) and (0 e - < 20° or 6 e - > 160°). 

The situation is completely different at .B-factories, where one can al- 
most completely reduce the background coming from e + e~ pairs ( FigEJ), by 
rejecting the events with at least one charged lepton, electron or positron, 
in the detector. 



4. Conclusions 

We have constructed the Monte Carlo generator EKHARA, which sim- 
ulates the reaction e + e~ — > 7r + 7r~e + e~ based on initial and final state emis- 
sion of a e + e~ pair from e + e _ — > ir + ir~ production diagram. A detailed 
study of this process, as a potential background for a(e + e~ — > ir + ir~) mea- 
surement via radiative return method, shows that it can become important, 
when the experimental precision will reach 1%, or better. 
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